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ABSTRACT 

We examine the origin and evolution of the mass-metallicity relationship (MZR, M*-Z) for galaxies 
using high resolution cosmological SPH + N-Body simulations that include a physically motivated 
description of supernovae feedback and subsequent metal enrichment. We discriminate between two 
sources that may contribute to the origin of the MZR: 1) metal and baryon loss due to gas outflow, 
or 2) inefficient star formation at the lowest galaxy masses. Our simulated galaxies reproduce the 
observed MZR in shape and normalization both at z=0 and z=2. We find that baryon loss occurs due 
to UV heating before star formation turns on in galaxies with M& ar < 10 s M Q , but that some gas loss 
due to supernovae induced winds is required to subsequently reproduce the low effective chemical yield 
observed in low mass galaxies. Despite this, we show that low star formation efficiencies, regulated 
by supernovae feedback, are primarily responsible for the lower metallicities of low mass galaxies and 
the overall M*-Z trend. We find that the shape of the MZR is relatively constant with redshift, but 
that its normalization increases with time. Simulations with no energy feedback from supernovae 
overproduce metals at low galaxy masses by rapidly transforming a large fraction of their gas into 
stars. Despite the fact that our low mass galaxies have lost a majority of their baryons, they are still 
the most gas rich objects in our simulations due to their low star formation efficiencies. 
Subject headings: galaxies: evolution — galaxies: formation — methods: N-Body simulations 



1. INTRODUCTION 

The observed trend that the metallicity of nearby 
galaxies decreases with galaxy m ass has been well es- 
tablished in th e literature (e.g., iTremonti et al.1 [2004t 
iLee et all 120061 ) . This tr end has also recently been ob- 
serve d at high redshifts (jSavag lioet al.l [20051 : lErb et all 
120061 hereafter S05 and E06, respectively). The gas 
metallicity of a galaxy depends on reprocessed material 
from its stars and subsequent inflow and outflow of gas. 
The mass-metallicity relationship (MZR) thus provides 
insight into two physical mechanisms important to the 
evolution of galaxies: star formation efficiency and gas 
inflow / outflow. 

In a "closed box" system (no inflow or outflow of ma- 
terial), metallicity Z is given by 



Z 



gas' 



(1) 



where y is the nucleosynthetic metal yield produced by 
stars and f gas is the gas mas s fraction (M ffas /[M gas 
+ M s tars]) (e.g.. iTinslevi Il980h . In general, low mass 
galaxies are observed to have higher gas f r actions tha n 
high mass galaxies (e.g., iGeha et alJl2006t IWestJ l2005h . 
This observation can easily explain the origin of the 
M*-Z relation (MZR) in the context of the closed 
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box model. In this scenario, low mass galaxies are 
gas rich because they are inefficient at turning gas 
into stars, likely due to their low surface densities 
(lKennicuttJll998HMartin fc Kennicu"ttll200H rVerde et all 



2002t iDalcanton et al.ll2004fl . 

Alternatively, SNe feedback might expel gas in galactic 
winds. In this scenario the loss of metals becomes pro- 
gressively more important at lower mass es due to shal- 
lower halo potential wells (for a review see lVeilleux et all 
l2005h and leads to the observed M*-Z trend. However, 
the importance of supernova feedback as a function of 
galaxy mass has been investigated by various groups , 
with results that are still debated (jVeilleux et alJ l2005t 
IGeha et alJl2006t lDalcantonll2006t E06). 

Recently, studies have investigate d how the MZR 
evolves with cosmic time (S05. E06. jDave et ail 120061: 
De Rossi et al.ll2006t iKobavashi et al.ll2006l : iTassis et al.1 



2006). Fully cosmological numerical experiments are one 



of the best tools to infer the relative importance of star 
formation (SF) and feedback processes in a realistic set- 
ting (e.g.. lOkamoto et "all 120051: iRobertson et all 120061 : 
iGovernato et alJ I2006L hereafter G06). In this Letter, 
we investigate the origin and evolution of the MZR in 
a concordance ACDM cosmological setting with N-body 
+ SPH simulations which include a physically motivated 
treatment of SF and SN feedback. These simulations al- 
low us to study the formation history of galaxies over 
range of masses, and thus quantify the importance of 
feedback and blowout versus low SF efficiencies over their 
lifetimes. 

2. THE SIMULATIONS 

The simulations used in this study are the culmina- 
tion of an effort to create realistic disk galaxies. As 
discussed in G06, at z=0 disk galaxies in our Simula- 
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tions fall on the Tully-Fi s her (TF) and baryoni c TF 
relations (jGiovanelli et al.l 119971: iMcGaughl I2005T) and 
the age- mass relation (|MacArthur et al.l 12004^ This is 
and important improvement over previous works, and 
is made possible due to the inclusion of a simple but 
physi cally motivated recipe (full details in lStinson et al.l 
2006) to describe SF and the effects of subsequent 
SNe feedback. O ur adopted SN feedbac k and cosmic 
UV background (|Haardt k. Madaul I1996D also drasti- 
cally reduce the number of galaxy satellites containing 
a signific ant stellar popula t ion, making many of them 
"dark" (jQuinn et alJ 119961 : iMoore et alJ 1 19991 ). Thus, 
the SN and SF efficiency parameters were adop ted as 
in G06, with an IMF from iKroupa et all (|1993[) . with 
no additional tuning of the parameters to affect met- 
als. Metal enrichment from both SN Iafcl l is fo llowed 
based on the pr escription of iRaiteri et all (|1996f ). who 
adopt yields fromlThielemann et al.l (|1986l) for SN la, and 
IWeaver fc Wooslevl ([19931) for SNII. 

We selected four "field" regions (dp/ p ~ 0.1) from a low 
resolut ion, dark matter (DM ) simulation run in GASO- 
LINE (|Wadslev et al.l [20041 using a concordance, flat, 
A-dominated cosmology: ilo = 0.3, A=0.7, h = 0.7, 
an = 0.9, shape p aram eter T = 0.21, an d f2& = 0.039 
(|Perlmutter et all 119971 : lEfstathiou et al.l l2002f) . These 
field regions were resimulated at higher r esolution using 
the v olume renormalization technique (Katz & White 
1993), which allows us to resolve fine structure while 
capturing the effect of large scale torques. The field re- 
gions are centered around galaxies with peak rotation 
velocity V ro t ~ 150 and 275 km/s, with each high res- 
olution region having a comoving volume of ^100-1000 
Mpc 3 . Two of the main galaxies are described in de- 
tail in G06, but for this study they have been rerun at 
higher resolution. At z=0, these runs yield a sample of 
31 high resolution disk and dwarf galaxies over a range 
of peak velocities 30 km/s < V ro t < 275 km/s, or total 
halo masses from 3.4 xlO 9 M Q to 1.1 xlO 12 M Q . Seven of 
these 31 galaxies are satellites of the primaries. Galax- 
ies and their parent halos were id entified using AHF 8 
(|Knebe et al.ll200lt IGill et al.ll200l . Our highest (low- 
est) resolution runs have particle masses of 7.6xl0 5 M Q 
(4.9xl0 6 M Q ), 1.3xl0 5 M (8.5xl0 5 M ), and 3.8xl0 4 
M (2.5xl0 5 Mq) for DM, gas, and stars, respectively, 
and a force resolution of 0.15-0.3 kpc. 

3. THE ORIGIN OF THE MZR 

The MZR for 31 high resolution galaxies in our simula- 
tions is shown in Fig.[TJ at both z=0 and z=2. The solid 
curved line is the empirica l fit fo r > 53,000 galaxies in 
SDSS from lTremonti et alJ (|2004 hereafter T04). It has 
been shifted down by 0.26 dex as found by E06 to cor- 
rect for the fact that the method used by T04 is known 
to systematically find O abundanc es up to 0.3 dex larger 
than other methods (T04, S05, E06 lLee et aT1l2006f ). The 
galaxies in our sample at z=0 match the empirical fit of 
T04 extremely well, reproducing the observed trend over 
almost 4 orders of magnitude in stellar mass. We note 
that our lowest mass galaxies, however, have a lower O 
ab undance than th at found for dwarf irregular galaxies 
bv lLee et al.l (|2006l) . 

8 AMIGA'S Halo Finder, av ailable for download at 
htt p : / / www .aip.de/P eople / aknebe /AMIGA] 



Fig.[T]also shows that the MZR for our simulated galax- 
ies at z=2 (open diamonds) is in excellent agreement with 
high redshift data. The error bars overplotted on Fig. [1] 
display the observational z=2 results from E06. We have 
identified in the simulations the z=0 galaxies correspond- 
ing to those plotted at z=2. The dotted lines in Fig. [1] 
connect a few of these galaxies, which are representative 
of how all the galaxies have evolved in the M*-Z plane 
with time. 

To make a meaningful comparison with the T04 sam- 
ple, Fig. [1] uses only cool gas particles (with temperature 
< 40,000 K) within the radius of each galaxy that en- 
closes 2/3 of the stellar mass. The median "observed" 
diameter for our galaxies is 5 kpc, similar to that found 
by T04. At z=2 we included all cold gas in each galaxy. 
Even if we include all cool gas at z=0, the M*-Z trend 
is unchanged up to 10 10 M Q , showing that the evolution 
seen in Fig. [T] is not an effect of our measurements. 

A few runs were repeated at lower resolutions (1/2 and 
1/20 of the mass resolution). Low resolution runs un- 
dergo less SF, thus decreasing both M* and the total 
amount of O produced. Our results from each run be- 
gin to significantly diverge from each other (AM* > 20%, 
Al2+log(0/H) > 0.2 dex) if the number of DM particles 
in a halo drops below a few thousand and/or mass and 
force resolution degrade too much. Note that this low 
resolution effect can lead to an artificially steeper M*- 
Z trend, making it essential that results converge. In- 
formed by these resolution tests, we have only included 
galaxies with more than 3500 DM particles in our sample 
(at high and low redshift), ensuring that their stellar and 
metal content have converged. 

Having shown that our simulations reproduce the ob- 
served MZR at both z=0 and z=2, we will next examine 
the origin of this relation. We focus on the importance 
of preferential metal loss versus low SF efficiency at low 
masses in producing the MZR. 

3.1. The Role of Mass Loss 

An observable quantity commonly used to investigate 
mass loss from galaxies is the effective yield, yeff- From 
Eq. [TJ we can define the effective yield, 

Veff = Z/(lnf- a l). (2) 

In a closed box system, y e ff is equal to the true stellar 
nuclcosynthctic yield, but it can be inferred from Eq. [2] 
that y e ff may deviate from the true yield due to inflow 
or outflow of gas. The top panel of Fig. [5] shows y e ff for 
our galaxies, derived again using inner cold gas abun- 
dances, versus the total baryonic mass of our galaxies at 
z=0 (black circles). The effective yields reproduce the 
observed trend from T04 (solid line) over most of our 
mass range. Our t wo lowest mass g alaxies have y e ff be- 
low that found bv lLee et al.l (|2006[ ). likely because these 
two have undergone additional tidal stripping of gas. 

The black circles in the bottom panel of Fig. [5]show the 
baryonic to dark mass ratio (B/DM) of our galaxies, nor- 
malized by the cosmic ratio (the seven satellite galaxies 
have been excluded from this panel). In our simulations, 
large deviations below B/DM = 1 can only be produced 
by mass loss. Our lowest mass galaxies have lost up to 
98% of their baryonic mass. 

Using the halos defined by AHF (our halo finder), we 
identify a given halo and its gas particles back to z=3. 
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Once all the gas particles that have ever belonged to a 
galaxy back to z—3 are identified, those particles that are 
no longer within the virial radius at z=0 are considered 
to be gas "lost" from the galaxy. The grey points in 
both panels of Fig. [2] represent this information, for all 
gas (remaining + lost). 

When we include all of the gas lost since z = 3, we 
find that our low mass galaxies still have B /DM that fall 
significantly below the cosmic ratio. The bottom panel 
of Fig.[2]show that our lowest mass galaxies have lost 50 - 
90% of their baryons prior to z = 3. We verify that a UV 
only simulation (SNc feedback turned off) removes up to 
80% of baryons from our lowest mass galaxies, but that 
these galaxies evolve as closed boxes once SF commences. 
We thus find that SNe induced mass loss can account 
for up to 20% of baryon loss from galaxies with M-bar 
< 1x10 s M©. Additional gas loss can occur via tidal 
stripping. 

Thus, Fig. [2] clearly demonstrates that gas has been 
preferentially lost from our lowest mass galaxies, and 
that this gas loss is responsible for the observed y e ff- 
Mf,ar trend. However, when all gas (remaining + lost) 
is used to rederive the MZR (Fig. [T]) , the observed trend 
that lower mass galaxies have lower metallicities still ex- 
ists, and it cannot be mass loss alone that causes the 
MZR. 

3.2. Star Formation Efficiencies 

As explained above, the mass loss experienced by our 
lowest mass galaxies cannot account for the overall M*- 
Z trend seen in Fig. [TJ We now investigate whether SF 
efficiencies can explain the origin of this trend. If lower 
mass galaxies convert gas into stars less efficiently than 
higher mass galaxies, then they will require a much longer 
period of time to enrich the remaining gas to the same 
metallicity as higher mass galaxies. This effect can lead 
to the observed MZR, with or without mass loss as an 
additional contributor. 

We verified that all of our ga laxies follow the 
Schmidt-Kennicutt law for SF (jKennicuttl 119981 : 
iMartin fc Kennicuttll200i[ ) . Our galaxies span a factor of 
250 in gas surface densities, yielding a SFR that varies 
by over three decades from our lowest to highest mass 
galaxies. Our lowest mass galaxies are thus incredibly 
inefficient at SF compared to our high mass galaxies. 
The consumption rate (SFR/M co /d gas ) of our lowest mass 
galaxies is about two decades lower than our highest mass 
galaxies. That is, while our highest mass galaxies will use 
up all their gas in 1 Gyr or less, our lowest mass galax- 
ies will take up to 60 Gyr if they continue SF at their 
current rate. 

Fig. [3] shows the gas mass fraction, f gas , versus ro- 
tational velocity for our galaxies (satellite galaxies ex- 
cluded) . As found by r ecent observation s of isolated 
galaxies (e.g. IWestl[2005l iGeha et al.ll2006[ h lower mass 
galaxies have increasing gas fractions, implying that 
lower mass galaxies have been less efficient at making 
stars. Note that while the lowest mass galaxies in our 
sample have lost more than 90% of their baryons (Fig. [2|) , 
they remain gas rich due to their low SF efficiencies (with 
the exception of two that have been tidally stripped of 
much of their gas) . 

We reran our simulations without SNe feedback. In 
our models, without this feedback the MZR is flat be- 



cause low mass galaxies experience unrealistically high 
SF rates, resulting in an overproduction of metals. This 
difference with the runs that do include SNe demon- 
strates that energy injection by SNe into the ISM reg- 
ulates SF efficiency, independent of whether SNe also 
lead to mass loss. SNe feedback thus indirectly affects 
the shape of the MZR by modulating the gas surface 
densities, and hence the SFR, with galaxy mass. 

4. THE EVOLUTION OF THE MZR 

Fig.[3]shows that our simulated galaxies at z=2 are less 
evolved than their z=0 counterparts, with z=2 galaxies 
being much more gas rich at a given rotational velocity 
than at z=0. By z=0, they have consumed more of their 
gas, increasing both their stellar mass and O abundance; 
hence the evolution of the MZR as shown in Fig. [1] 

Both S05 and E06 examined the MZR at z~0.7 and 
z^2, respectively, and drew conclusi o ns sim ilar to ours. 
As pointed out by S05 and iPettinil (|2006f) . at a given 
metallicity the high z data are about an order of mag- 
nitude more massive than the z=0 trend (Fig. [J). This 
suggests that more massive galaxies enrich at a faster 
rate than lower mass galaxies. This is supported by our 
results in section [3~2l Both S05 and E06 conclude that 
the evolution of the MZR is best described by a model in 
which low mass galaxies are inefficient at SF relative to 
high mass galaxies for a similar formation epoch. Thus 
the period of SF is lengthened in lower mass galaxies and 
they remain more gas rich at a given redshift. 

We find that a shift downward by 0.3 dex in Alog(0/H) 
from fit shown in Fig.[T]is an excellent fit to our z=2 data, 
as was also found by E06. A shift in Alog(M») is much 
too steep to match our simulations. Thus, we find that 
the slope of the MZR remains relatively constant with 
redshift, but the normalization increases with time. 

5. CONCLUSIONS 

We have presented results from high resolution simu- 
lations of field galaxies in a ACDM context. The sim- 
ulations are carried to z=0 and include SF and a phys- 
ically motivated description of the effects of SNe feed- 
back and metal enrichment. We have shown in previous 
work (G06) that this approach leads to disk galaxies that 
fall on the TF and baryonic TF relations and reproduce 
the age-mass relation of disk galaxies and abundance of 
galaxy satellites. Here we present a natural extension 
of this analysis, showing that the simulated galaxies are 
also in excellent agreement with the observed MZR both 
at z=0 and z=2. 

Simulated galaxies match both the slope and normal- 
ization of the MZR over four orders of magnitude in 
stellar mass. We find that the slope of the MZR is 
roughly constant with redshift, but that its normaliza- 
tion increases with time. Such evolution agrees with the 
high-z obse rvations of E06, and with early theo retical 
predictions (|Dave et al.ll2006l : iDe Rossi et alj|2006f) . Fur- 
thermore, the mass scales of the turnovers in our sim- 
ulated Z, y e ff, B/DM, and f gas versus mass relations 
are in good agreement with a range of previous observa- 
tional and theoretical studies that find critical transitions 
at ~120 km/s (e.g., iDalcanton et al.l l 2004t iKeres et all 
[20051 iDekel fc Birnboim 2006). 

Galaxies with M bar < 10 8 M© (V rot ~ 50 km/s), lose 
a significant amount of their baryons due to combined 
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heating from the cosmic UV held and SN feedback. 
Baryon loss due to SN feedback results in a y e ff-Mb ar 
trend in agreement with the results from the SDSS (T04). 
However, including gas that has been lost from galax- 
ies does not alter the O abundances significantly, and 
does not change the overall shape of the MZR, confirm- 
ing that the lower O abundances in our low mass galaxies 
are primarily due to their low SF efficiencies rather than 
directly to blowouts. Despite strong baryon loss from 
their halos, we find that our cold gas fractions at low 
galaxy masses are in qualitative agreem ent with recent 
observational samples (|Geha et~a l. 2006). 

SN feedback plays a crucial role in lowering the SF effi- 
ciency in low mass galaxies and originating the turnover 
in the MZR at the observed scales. Without energy in- 
jection from SNe to regulate SF, gas that remains in 
galaxies rapidly cools, forms stars, and increases the O 
abundance, even by z=2 (see G06 for additional short- 
comings of the no feedback case). With feedback from 
SNe turned off, small galaxies produce too many metals 
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Fig. 1. — The mass-mctallicity relation of our simulated galaxies at z=0 (solid points) and z=2 (open diamonds). The solid curved line is 
the observational fit to >53,000 galaxies in SDSS from T04, shifted down by -0.26 dcx as found by E06. Error bars show the observational 
mass-metallicity relation at z=2 from E06. Dotted lines connect some of the z=0 galaxies to their progenitors at z=2, showing how galaxies 
evolve in the M*-Z plane with time. 



6 



Brooks et al. 





— 1.3 




-2.0 


'a; 


-2.5 


cn 




o 
_i 


-3.0 




-3.5 




-4 




1.0 






Q 




m 


0.10 



0.0 



10 




10" 10" 10 

Baryonic Mass (M Q ) 



Fig. 2. — Top Panel: Effective yield vs total baryonic mass for our galaxies at z=0 (black circles). The solid curved line is the 
observational fit to >53,000 galaxies from T04. Bottom Panel: Baryonic to dark matter mass ratio (relative to the cosmic abundance) vs 
total baryonic mass. Satellite galaxies excluded. In both panels, grey circles show what the value would be if no gas loss from the galaxies 
had occurred. 
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Fig. 3. — The gas mass fraction vs Vrot for our galaxies at z=0 and z=2 (satellite galaxies excluded). Faint dotted lines connect a few 
of the galaxies at z=2 to their evolved counterpart at z=0. 



